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Abstract 

For two or more classes of points in R'^ with d > 1, the class cover catch digraphs (CCCDs) can 
be constructed using the relative positions of the points from one cla ss with respect to the points from 
the other class. The CCCDs were introduced bv lPriebe et alj l|200ll ') who investigated the case of two 
classes, X and y . They calculated the exact (finite sample) distribution of the domination number of the 
CCCDs based on X points relative to y points both of which were uniformly distributed on a bounded 
interval. We investigate the distribution of the domination number of the CCCDs based on data from 
non-uniform X points on an interval with end points from y . Then we extend these calculations for 
multiple y points on bounded intervals. 
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1 Introduction 



In 2001, a new classificatio n method was deve loped which was based on the relative positions of the data 
points from various classes: iPriebe et al.l (j200l[ ) introduced the class cover catch digraphs (CCCDs) in K and 
gave the exact distribution of the domin ation number of the CCCDs for two classes, X and 3^, with uniform 
distribution on a bounded inter val in R. DeVinnev and Wierman ( 2003) proved a SLLN result for the one - 
dim ensional class cover pr oblem. DeVinnev et al.l ( 2002 ). Marchette and Prieb3 ( 2003f ). Priebe et al. ( 2003a ). 



and lPriebe eTal] (|2003b[ ) extended the CCCDs to higher dimensions and demonstrated that CCCDs are a 
competitive alternative to the existing methods in classification. The classification method based on CCCDs 
involves data reduction (condensing) by using approximate — rather than exact — minimum dominating 
sets as prototype sets, since finding the exact minimum dominating set for CCCDs is an NP-hard problem 
in general. However for fi n ding a dominating set of CCCDs on the real line, a simple linear time algorithm 
is available ( Priebe et al.l ( 200lf l). But unfortunately, the exact and the asymptotic distributions of the 
domination number of the CCCDs are not analytically tractable in multiple dimensions. 



To addr ess the latter issue of intractability of the distribution of the domination number in multiple 
dimensions, Cevhan and Priebd ([2003., ,20051 ) introduced the central similarity proximity maps and r-factor 
proportional-edge proximity maps and the associated random proximity catch digraphs. Proximity catch 
digraphs are a generalization of the CCCDs. The asymptotic distribution of the dominati on number of the lat- 
ter is c alculated and then used in testing spatial patterns between two or more classes. See lCeyhan and Priebe 
( 20051 ) for more detail. 



In this article, wc generalize the original result of lPriebe et al.l (j200l[ ) to the case of non-uniform X points 
with support being the interval with end points from 3^, and then to multiple y points in a bounded interval 
(c, d) C M with c < d. These generalizations will also serve as the bases for extension of the results for the 
uniform and non-uniform data in higher dimensions. 



2 Data-random Class Cover Catch Digraphs 



Let be a measurable space and Xn — {Xi,...,X„} and J^m = {yi: ■ ■ ■ ,ym} be two sets of re- 

valued random variables from classes X and y, respectively, with joint probability distribution Fx.y- Let 
•) : X [0, oo) be any distance function. The class cover problem for a target class, say X , refers to 
finding a collection of neighborhoods, Ni around Xi such that (i) Xn C {UiNi) and (ii) J^m H {y}iNi) = 0. 
A collection of neighborhoods satisfying both conditions is called a class cover. A cover satisfying condition 
(i) is a proper cover of class X while a cover satisfying condition (ii) is a pure cover relative to class 3^. This 
article is on the minimum cardinality c lass covers; t h at is, class covers satisfying both (i) and (ii) with the 
smallest number of neighborhoods. See lPriebe" et all (120011 ). 



Consider the map N : Vl 2^^ where 2^ represents the power set of Vt. Then given y,n Q f^, the 
proximity map Ny{-) : f2 — *■ 2^ associates with each point x G fl a. proximity region Ny{x) C 57. For _B C 17, 
the Fi-region is the image of the map ri(-, Ny) : 2^ ^ 2^ that associates the region Ti^B, Ny) := {z e i7 : 
B C Ny{z)} with the set B. For a point a; G 17, we denote Fi({x}, Ny) as Fi(a;, Ny). Notice that while the 
proximity regions are defined for one point, Fi-rcgions are defined for sets of points. 



The data-random CCCD has the vertex set V = Xn and arc set A defined by {Xi,Xj) E A 
Xj £ Ny{Xi). In particular, wc use Ny{Xi) = B(Xi,ri), the ope n ball around Xi with radius 



TmTiY^y^d{Xi,Y), as the proximity map as in iPriebe et al.l (|200l[ ). We call such a digraph a 
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digraph. A ^n m-digraph is a pseudo digraph according some authors if loops arc allowed 
Chartrand and Lesniak 



A data-random CCCD for il = M'' and A^^ = B{Xi,ri) is referred to as "^n^m-graph in Priebe et al. (12001 ). 
We ch ange the notation to emphasize the fact that S>n,m is a digraph. Furthermore, ICevhan and Priebe 
(2003) call the proximity map A^^ = B{Xi, ri) a spherical proximity map. 



The ^„_m-digraphs are closely related to the proxi mity graphs o f I Jaromczvk and Toussaintl (Il992l) and 



might be considered as a special case of covering sets of iTuzal (|1994I ) and intersection digraphs of ISen et al 



( 1989|). Our data-random proximity digraph is a vertex-random proximity digraph and not a standard one (see 



e.g., Janson et al. The randomness of a ^„^m-digraph lies in the fact that the vertices are random 

with the joint distribution Fx,y, but arcs {Xi,Xj) are deterministic functions of the random variable Xj 
and the random set A^,-. 



3 Domination Number of Random ^„ ^-digraphs 



In a digraph D = (V, A) of order | V| = n, a vertex v dominates itself and all vertices of the form {u : (v, u) € 
A}. A dominating set, Sd, for the digraph _D is a subset of V such that each vertex v G V is dominated 
by a vertex in So- A minimum dominating set, S}), is a dominating set of minimum cardinality; and the 
domination number, denoted j{D), is defined as j{D) := \S'^\, where | ■ | is the set cardinality functional 



I WestI ((20011)). If a minimum dominating set consists of only one vertex, we call that vertex a dominating 
vertex. The vertex set V itself is always a dominating set, so j{D) < n. 



I'' with P{X 

and ym arc random samples from Fx and Fy, respectively. 



Let {W') := {Fxx on 
^„,m-digraphs for which X, 

distribution of X, Y is Fx,y € (M^) . We call such digraphs as T 

on the random variable "f{D). To make the dependence on sample sizes explicit, we use "f{D. 
7(D). It is trivial to see that 1 < ^{Dn,rn) < n, and j{Dn^m) < n for nontrivial digraphs. 



y) = 0}. As in lPriebe et al.l (j200lh . in this article, we consider 

and the joint 
random Sln.m- digraphs and focus 
) instead of 



4 The Distribution of the Domination Number of 
^„„,-digraphs 



-random 



In 



th e dat a-random CCCD is a special case of interval catch digraphs (see, e.g.. Sen et al. ( 1989f ) and 
Prisnei ( 1994 )). Let Xn and be two samples from and Y(j) be the order statistic of for 

j ~ 1,2, . . . ,m. Then Y(^j) partition M into (yn + 1) intervals. Let 



< ^(1) < • • • < ^(m) < ^(m+l) := 00, 



and X, 



U)) 



X^ 



Xn n Ij, and := {y(j_i), F^)} for j = 1,2,..., (m + 1). This yields a 
disconnected digraph with subdigraphs for j ~ 1,2,..., (m + 1), each of which might be null or itself 
disconnected. Let ^{D^) denote the the cardinality of the minimum dominating set for the component of 
the random ^„^m-digraph induced by the pair X^ and y^ , nj \X^, and Fj be the density Fx restricted 
to Jj. Then 7(£'„,,n) = T^. J^i^ l(F>^). ^ye stu dy the simpler random variable j{D^) first. The following 
lemma follows trivially fsee fPriebe et al.l (|200l[ )). 

Lemma 4.1. For j g {1, (m -I- 1)}, ^{D^) = l{nj > 0) where I(-) is the indicator function. 
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For j = 2, . . . , m and rij > 0, we prove that ^{D^) G {1, 2} with the distribution dependent probabiUties 
^~Pn {Fj),Pn {Fj), respectively, where — P{'y{D^) ~ 2). A quick investigation shows that 7(D-') = 2 

iff X^n (fllf!ik)±!(2^, min(A'^)+yo) ^ J ^, ^j^^^ .^^ ^ B{x,r{x)) iff Z S ( "'^^^ (^^H^,-!) ^ mhU£i)+lo,^ 

where r(x) = min(a; - Y(j_r),Y(j) - x). Hence Ti{XJ,Ny) = |^ '"ax(^^)+y(,_i) ^ min(^^)+yo) ^ ^ 
definition, if A'-' n Vi{X^ , Ny) ^ 0, then j{D^) = 1; hence the name Ti-region and the notation ri(-, iV^;). 
Theorem 4.2. i^or j — 2, . . . ,m, jiD^) ^ 1 + Bernoulli (Fj)) for rij > 0. 



Proof: See iPriebe et alj (|200l[ ) for the proof 



The probability P{"f{D^) = 2) = P{X^ n Vi[X' ,Ny) = 0) depends on the conditional distribution Fx\y 
and the interval Ti[X^ , Ny)^ which, if known, will make possible the calculation of p„ [Fj). As an immediate 
result of Lemma I4?l1 and Theorem 14. 2[ we have the following upper bound for 7(-D„,m)- 

Theorem 4.3. Let -D^.m be an J-{W)-random ^n.m-digraph with n > 0, m > and ki and k2 be two natural 
numbers defined as ki I]Jl2 Id-^nnX, ! > 1) and k2 Z]j=2 Id-^nnXjl = l)+J2je{i,{m+i)} ^i^nf^^j ^ 0)- 
Then 1 < ^(Dn.ni) < 2 ki + k2 < min(7i, 2 m). 



In the special case of fixed 3^2 = {yi7 72} and Xn a random sample from W(yi, y2), the uniform distribution 
on (yi,y2), we have a ^„_2-digraph for which Fx = W(yi,y2) and Fy is a degenerate distribution. We call 
such digraphs as ZY(yi, y2)-ran(iom ^^n, 2 -digraphs and provide an exact result on the distribution of their 
domination number in the next section. 



4.1 The Exact Distribution of the Domination Number of W(yi, y2)-random ^„,2- 
digraphs 



Suppose 3^2 = {717X2} C M with —00 < yi < y2 < 00 and Xn = {Xi, . . . , Xn} a set of iid random 
variables from U{yi,y2)- Any Z^(yi,y2) random variable can be transformed into a U(0, 1) random variable 
by (j){x) = {x — yi)/(y2 — yi), which maps intervals (ti,t2) C (yi,y2) to intervals (0(<i), (I){t2)) ^ (0, 1). So, 
without loss of generality, we can assume Xn = {Xi, . . . , Xn} is a set of iid random variables from the U (0, 1) 
distribution. That is, the distribution of j{Dn,2) does not depend on the support interval (yi,y2)- Recall 
that -/(Dn.2) = 2 iff Xn n ri(A'„,7Vy) = 0, then P{-f{Dna) = 2) = 4/9 - (16/9)4-". For more detail, see 
(jPriebe et al.l (|200lh ). Hence, for U{yi,y2) data, wc have 



7(^n.2) 



1 w.p. 5/9+ (16/9) 4- 

2 w.p. 4/9- (16/9) 4- 



for all n > 1, 



(1) 



where w.p. stands for "with probability". Then the asymptotic distribution of ^{Dn.2) for hi{y 1,^2) data is 
given by 

1 w.p. 5/9, 

2 w.p. 4/9. 



lim 7(Ai,2) 



(2) 



For TO > 2, iPriebe et al.l (j200l[ ) computed the exact distribution of 7(-D„.m). However, independence of 
the distribution of the domination number from the support interval docs not hold in general; that is, for 

Xi '~ F with support S{F) C (yi,y2), the exact and asymptotic distribution of ^{Dn.2) will depend on F 
and 3^2 ■ 
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4.2 The Distribution of the Domination Number for jF(M)-random ^„ 2-digraphs 

For 3^2 = {yi,y2} C K with — oo < yi < y2 < oo, a quick investigation shows that the Fi-region is 
Ti^Xn, Ny) ~ ^ yi+^c) ^ y2+^£(i)_~j ^ ]\jote that n ri(A'„, Ny) is the set of all dominating vertices, which is 

empty when 7(-D„.2) > 1- To make the dependence on F explicit and for brevity of notation, we will denote 
the domination number of the -F((yi, y2))-random ^„_2-digraphs as 7„(i^). 

Let Pn{F) P{jn{F) = 2) and p{F) := lim„^oo P(7„(F) = 2). Then the exact (finite sample) and 
asymptotic distributions of jn{F) are 1 + Bernoulli (p„(F)) and 1 + Bernoulh {p{F)), respectively. That is, 
for finite n, we have 

7n{F)^l I "'■P- foraUn>l. (3) 

'"^ ^ \ 2 w.p. pniF) - ^ > 

The asymptotic distribution is similar. 

With 3^2 = {0, 1}, let be a distribution with support S(F^ C (0, 1) and density / and let Xn be a set 
of n iid random variables from F . Since ^n(F) £ {1, 2}, to find the distribution of -fn{F), it suffices to find 
P{^^{F) 1) or P{jn{F) = 2). For computational convenience, we employ the latter in our calculations. 

Then 

JS(F)\ri(X„,Ny) 

where /i„(xi,x„) = n{n — 1) [F{xn) — ^(a^i)]" fi^i) fi^n) 1(0 < < x„ < 1) which is the joint 
probability density function of ^(i), ^(n)- 

If the support S{F) ~ (0, 1), then the region of integration becomes 

|(a;i,a;„) e (0, 1)^ : (1 + a;i)/2 < x„ < 1; < xi < 1/3 or 2 < x„ < 1; 1/3 < zi < 1/2|. 

The integrand in Equation ([4|) simplifies to 

H{xi , Xn) := n (n - l)/(xi)/(a;„) [Fixn) + F (a;„/2) - {F ((1 + xi)/2) + F{xi))] (5) 

Let Xn be a set of iid random variables from a continuous distribution F with S{F) C (0, 1). The simplest 
of such distributions isZY(0, 1), the uniform distribution on (0, 1), which yields the simplest exact distribution 
for 7„(F). If X ~ F, then by probability integral transform, F{X) ~ W(0, 1). So for any continuous F, we 
can construct a proximity map depending on F for which the distribution of the domination number for the 
associated digraph will have the same distribution as that of 7„(W(0, 1)). 

Proposition 4.4. Let Xi F which is an (absolutely) continuous distribution with support S{F) = (0, 1) 
and Xn = {Xi,...,Xn}. Define the proximity map Nf{x) := F-'^{Ny{F{x))) = F-^{B{F{x),r{F{x)))) 
where r(F{x)) ~ min(i^(x), 1 — F{x)). Then the domination number of the digraph based on Np, Xn, and 
3^2 = {0, 1}, is equal in distribution to 7„(Z^(0, 1)). 

Proof: Let Ui := F{Xi) for i = l,...,n and Un = {Ui, . . . ,Un}- Hence, by probability integral 

transform, Ui Z^(0, 1). Let [/(fc) be the fc*'' order statistic of Un for fc = l,...,7i. Furthermore, such 
an F preserves order; that is, for x < y, F{x) < F{y). So the image of Nf{x) imdcr F is F{Np{x)) = 
Ny{F{x)) = B{F{x),r{F{x))) for (almost) aU x g (0,1). Then F(7Vf(X,)) Ny{F{X,)) = Ny{U,) for 



1 - 



F{{l+xi)/2)-F{x„/2) 



F{x„)-F{xi) 



fin{xi,x„)dx„dxi, 



(4) 
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i = 1, . . . ,n. Since Ui Z-/(0, 1), the distribution of the domination number of the digraph based on Ny, 
Un and {0, 1} is given in Equation (P). Observe that Xj e Nf{X,) iff G F-^{B{F{X,),r{F{X,)))) iff 
F{Xj) e B{F{X,),r{F{X,))) iff U, e B{U^,r{U^)) for ^,J = f,...,n. Hence P{Xn C A^F(^^)) = -P(Z^» C 
Ny{Ui)) for aU i = 1, . . . , n. Therefore, n ri(A'„, TVi^) = iff n Ti(Un, Ny) = 0, which imphes that the 
domination number of the digraph based on Np, Xn, and 3^2 = {0, f } is 2 with probabihty 4/9 — (16/9) 4~". 
Hence the desired resuh follows. ■ 



For example for F{x) = 1(0 < x < 1) + I(a; > 1), 



{0,V2x) forxe [0,1/V2] 

(V2a;2 - f,l) for x e (l/V2,l] 



There is also a stochastic ordering between jn{F) and 7„(Z//(0, 1)) provided that F satisfies some condi- 
tions which are given in the following proposition. 

Proposition 4.5. Suppose Xn = {^i, ■ . ■ ,Xn} is a random sample from a continuous distribution F with 
S{F) C (0, 1) and let be the j*'* order statistic of Xn for j = 1, . . . ,n. If 

F(X(„)/2) < F (X(„)) /2 and F{X^,^) < 2 F ((l + X^,^) /2) - 1 hold a.s., (6) 

then 7n(F) <'^-^ 7„(Z^(0, 1)). If <'s in expression ([6]) are replaced with >'s, then 7,i(F) >"^-^ 7„(Z^(0,1)). 

// < 's in expression ([6]) are replaced with ~ 's, i/ien jn{F) = 7„(Z^(0,1)) where = stands for equality in 
distribution. 



Proof: Let Ui :— F{Xi) for i = 1, . . . , n and Un = . ■ . , Un}- Then, by probability integral transform, 
Ui *~ U{0,1). Let J7q) be the j*'* order statistic of Un for j = l,...,n. The Fi-region for Un based 
on Ny is TiiUn,Ny) = ([/(„)/2, (l + ?7(i))/2); likewise, Ti{Xn,Ny) = (X(„)/2, (l + X(i))/2). But the 
conditions in expression ® imply that Ti{Un,Ny) C F(ri(A'„, A^j;)). So W„ n F(ri(A'„, A^j;)) = implies 
that Un n ri(W„, TVj;) = and W„ n F{Ti{Xn, Ny)) = iff A'„ n Ti{Xn,Ny) = 0. Hence 

Pn{F) = P(A'„ n ri{Xn,Ny) = 0) < P(Z^„ H Ti{Un, Ny) = 0) = p„(Z^(0, 1)). 

Then 7„(F) 7„(Zi(0, 1)) follows. The other cases can be shown similarly. ■ 

For m ore on the com parison of jn{F) for general F against 7„(Z//(0, 1)), see Section 4.2.2 of the technical 
report bv lCevhanI (|2004l ). 



4.2.1 The Exact Distribution of jn{F) for F with Piecewise Constant Density 

Let 3^2 ~ {0, 1}. We can find the exact distribution of jn{F) for F whose density is piecewise constant. Note 
that the simplest of such distributions is the uniform distribution U{0, 1). Below we give some examples for 
such densities. 

Example 4.6. Consider the distribution F with density /(•) which is of the form f{x) = y^jj i{S < x < 
1-6) with 6 e [0, 1/2). Then F{x) = l{8 < x < I - 8) + l{x > I - 5) . The integrand in Equation ^ 
becomes 

n(n-l) |^ 3(.T„-xi)-i y~^ 
(1-2 5)2 1^ 2(1-2(5) ) ■ 
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Then for S S [0,1/3] 

Pn{F) = / / H{xi,Xn)dXndxi + / / H {xi , Xn) dXndxi 

JS J{l+xi)/2 Jl/3 J2xi 

'1-35V 



= (4/9-(16/9)4-")(i-||) 



(7) 



which converges to as n ^ oo at (an exponential) rate 0((i^)"). For S e [1/3,1/2), it is easy to 
see that ^n{F) — 1 a.s. In fact, for 5 € [1/3,1/2) the corresponding digraph is a complete digraph of 
order n, since Xn C N{Xi) for each i = 1, . . . , n. Furthermore, if 5 = 0, then F = U{Q, 1) which yields 
p„(F) =4/9- (16/9)4-". □ 

Example 4.7. Consider the distribution F with density /(•) which is of the form 

fix) = I{x e (0, 1) \ (1/2 - S, 1/2 + S)) with S e [0, 1/6]. 

Then the cumulative distribution function (cdf) is given by 

F{x) = i^i(x)l(0 <x < 1/2-6) +F2 {x)l{l/2-d < X < 1/2 + 6) + F3 {x)I{l/2 + 6 < a; < l) +I(a: > l), 
where 

Fi{x) ^x/{l-2S), F2(x) = l/2, andF3(a;) = (x-2(5)/(l-2(5). 

There are four cases regarding the relative position of X(„)/2, (l + ^(1)) /2 and 1/2 — 6, 1/2 + S that yield 
IniF) = 2: 

case (1) (X(„)/2, (l + X(i))/2) C (1/2-5,1/2 + 5); case (2) X(„)/2 < 1/2 - 5 < (l + X(i))/2 < 1/2 + 5; 
case (3) 1/2 - 5 < X(„)/2 < 1/2 + 5 < (l + X(i))/2; case (4) X(„)/2 < 1/2 - 5 < 1/2 + 5 < (l + /2. 

Let Ej{n) be the event for which case (j) holds for j — 1, 2, 3, 4, for example, 

E^{7i) := {(X(„)/2,(l + X(i))/2) C (1/2-5,1/2 + 5)}. 

Thenp„(i^) = Y7j=iP {ln{F) 2, Ej{n)). Furthermore, cases (2) and (3) are symmetric; i.e., P(7„(F) — 
2,E2{n)) = P{ln{F) = 2, £'3(71)). Then in case (1), we obtain P{-f„{F) = 2,Ei{n)) = 1-2 (Mi)" + 

IZ2S ) ■ Note that P(ri(A'„, iV^;) C (1/2 — 5, 1/2 + 5)) ^ 1 as n ^ 00, hence it suffices to use this case 
to show that Pn{F) 1 as n cx) at an exponential rate since P{Ei{n)) < pn{F). 

In cases (2) and (3), we obtain P(7„(-F) = 2, E2{n)) = | (l - ^) _ (j^) ") ^nd in case 

(4), P(7„(F) = 2,Ei{n)) = | (1-4-"+^) (yEH)" • See lCevhanI (|2004l ) for the details of the computations. 

Combining the results from the cases, for 5 e [0, 1/6] we have 

P{^n{F) = 2) = 1 + (^\^y (1/9+ (32/9)4-) - ([^)" (2/3 + (16/3)4"") , (8) 

1-4(5 



which converges to 1 as n —> 00 at rate O I I j 



-2S 



Notice that if 5 = 0, then F — U{0,1). The exact distribution f or 5 g (1 / 6, 1/3 ) can be found in a similar 
fashion. Furthermore, if 5 G [1/3, 1/2], then Pn{F) = 1-25". See lCevhanI (|2004f l also for the details of the 
computations. □ 
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Example 4.8. Consider the distribution F with density /(•) which is of the form f{x) = (1 + S)l(^x G 
(0, 1/2)) + (1 - (5) I{x e [1/2, 1)) with S e [-1, 1]. 



Then 



See CevhanI ( 2004 ) for the derivation. Hence \min—,aoPn{F) = '^^g_g2 ^ —'■ Ppi^), with the rate of convergence 



O ^("'^)"^- Note that pf{5) G [0,4/9] is continuous in 5 and decreases as increases. If (5 = 0, then 
F = U{Q, 1) and pf{S = 0) = 4/9. Note also that pf{S = ±1) = 0. □ 

Example 4.9. Consider the distribution F with density /(•) which is of the form 

fix) = (1 + 5)1(0 < a; < 1/4) + (1-5)1(1/4 <x < 3/4) + (1 + 5)1(3/4 < a; < 1) with 6 G [-1,1]. 



The exact value of p„(i^) is available, but it is rather a lengthy expression (see ICevhanI (|2004[ ) for the 



expression and its derivation). But the limit is as follows: Pn{F) — > ^(3^^*) I Pf{5) as n — > 00 with the rate 
of convergence O (^(^^)"^- So pf{5) is increasing in 5. Notice here that Pn{F) and pf{S) are continuous 
in 5 and pf{S) > for all 5 e (-1, 1]. Moreover, pf{S = 1) = 1 and pf{S = -1) = 0. □ 

Note that extra care should be taken if the points of discontinuity in the above examples are different 
from {1/4, 3/4} or 1/2, since the symmetry in the probability calculations no longer exists in such cases. 



4.2.2 The Exact Distribution of jn{F) for Polynomial / Using Multinomial Expansions 

The exact distribution of ^n{F) for (piecewise) polynomial f{x) with at least one piece is of degree 1 or 
higher can be obtained using the multinomial expansion of the term (^"^^ hi Equation ([5|) with careful 
bookkeeping. However, the resulting expression for Pn{F) is extremely lengthy and not that informative. 

The simplest example is with f{x) ^2x and F{x) = x"^. Then Pn{F) = P{jn{F) = 2) = Ai{n) + A2{n), 

1/31 I/2X 

where Ai(n) := Ji--^^^^y^H{xi,Xn)dxndxi, A2{n) := J^^^ J^^^ H{xi,Xn)dxndxi, and H{xi,Xn) = 
nin — l)xi a;„ (5 — 1 — 2 xi — 5 a^i)" ^. Then 
/•1/3 

Klin) = / (8na;i/5) (l - a;i/2 - 5 .t2/4)"^^ - (Snxi/S) (I/I6 + a;i/2 - 15x^/16)""^ dcci. 
Using the multinomial expansion of (0"^^ with respect to xi in the integral above, we have 

. . ^ V^/ n-1 \ 8ri (-1)92+915-1+912-92-2913-2-92-2,1 j^(_2^U+9l 23-392-493-491 ;L5gi3-2-92-29i 

Al(71,) = > ^ \ ^ ^ 

j^\qi.q2.qz 2 + q2 + 2qi 5(2 + g2 + 2gi) 

where Q2 = {gi, 92, 93 e N : gi + 92 + 93 = n- l}. 

Similarly, the second piece follows as 
/■I 

A2 (n) = / {8nxi/5) (l - a:i/2 - 5 a;?/4)""^ - (Snxi/S) (l5/a;?/4 - 1/4 - xi/2)""^ dxi. 
^1/3 
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Again, using the multinomial expansion of the (•)" ^ term above, we get 
A2(n)=y f "'^^ V2nf9(-l)''^+''i5''M-2'-i-'-^ +9(-l)i+'''^+''n5''U 

90 + 45 r2 + 90 n 



^11. jy—j-j u 1 ^ ^_^'il+'r3+r'2 1 c;!-!/! -2 ri-r3-r2 



where Q3 = {ri, r2, ra e N : ri + 7-2 + = n — l}. See ICevharJ (|2004[ ) for more detail and examples 



For fixed numeric ri, one can obtain Pn{F) for F (omitted for the sake of brevity) with the above densities 
by numerical integration of the below expression. 

p„(F) = P(7„(F) = 2) = / / H{xi,Xn)+ / H{xi,Xn)dxndxi, 

Jo J{l+xi)/2 Jl/3 J2xi 

where H{xi, Xn) is given in Equation 

Recall the JF(R'')-random ^„.,„-digraphs. We call the digraph which obtains in the special case of 
y-m, = {yi!y2} and support of Fx in (yi,y2)7 J-{(yi,y2))-™ndom ^^n. 2 -digraph. Below, we provide asymptotic 
results pertaining to the distribution of such digraphs. 



5 The Asymptotic Distribution of the Domination Number of 
jr((yi, y2))-random ^„ 2-digraphs 

Although the exact distribution of %-i{F) is not analytically available in a simple closed form for F whose 
density is not piecewise constant, the asymptotic distribution of "fniF) is available for larger families of 
distributions. First, we present the asymptotic distribution of 7„(F) for ^„_2-digraphs with 3^2 = {yi7y2} C 
R with yi < y2 for various F with support S{F) C (yi,y2). Then we will extend this to the case with 
yrn C K for m > 2. 

For e G (0, (yi + y2)/2), consider the family of distributions given by 
^((yi,y2),e) = {f: (yi,yi+£)U(y2-£,y2)u((yi+y2)/2-£,(yi+y2)/2 + e) C 5(F) C (yi, y2)}. 



Let the fc*'* order right (directed) derivative at x be defined as /'■'^■'(x+) := lim;)_,o+ ~ ^'^'^^l ^ ^ 

for all fc > 1 and the right limit at c be defined as /(c+) := lim;i^Q+ /(c + K). The left derivatives 
and limits are defined similarly with +'s being replaced by — 's. Furthermore, let h = ijrix^h'i) and c. ~ 
(ci,C2) and the directional limit at (ci,C2) G M'^ for g{x,y) in the first quadrant in be g(cj'",cj) := 
lim||^ll^g g{c+ h) and the directional partial derivatives at (ci,C2) along paths in the first quadrant be 

/ll,/i2>0 

d''+^g{c+,c+) 1 (d^gjc+h) d^glc)' 



hi,h2>0 

Theorem 5.1. Let 3^2 = {yi,y2} C K with —00 < yi < y2 < 00 and Xn = {A^i,---,^ri} with Xi 
F € JF((yi, y2), e) . Let Dn,2 he the random &n.2-digraph based on Xn and 3^2- Suppose k > is the 
smallest integer for which F{-) has continuous right derivatives up to order [k + 1) at yi, (yi + y2)/2, 
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/('^')(y+) + 2-('=+i) Z^'^) ((^4^)"^) 7^ and f^^Hy^) = for all j = 0, 1, . . . , fc - 1; and £ > is the smallest 
integer for which F{-) has continuous left derivatives up to order {£ + 1) at y2, (yi + y2)/2, f^^\y2) + 
2-(c+i) |^(yi+y2)-^ ^ q and f^^\y2) = V a^^? = 0, 1, 1. T/ien 7„(F) - 1 + Bernoulli 

where Pn{F) :~ P{pin{F) = 2) and for hounded /^'^■'(•) and f^^\), we have the following limit 



lim PniF) 



/W(y+) + 2-(fc+i) /(fc) ((yi±y2)+)] [/^(y^-) + 2-(^+i) /W ((^) ) 



Note also that pi{F) = 0. 

Proof: First suppose (yi,y2) = (0,1). Recall that ri{Xn,Ny) = (X(„)/2, (l + X(i))/2) C (0,1) and 
7„(i^) = 2 iff A'„ n ri(A'„, Ny) = 0. Then for finite n, 



Pn{F)=P{pin{F)=2)= j H{xt,Xn)dXr 

Js{F)\ri{X„,Ny) 

where H{xi, Xn) is as in Equation ([5]). 



dxi, 



Let e € (0, 1/3). Then P(X(i) < e, >1 — e)^lasn— >oo with the rate of convergence depending 

on F. So for sufficiently large n, 

nl r in-2 

n{n~l)fixi)f{xr,) F{x,-,)~F{xi) + F{x,j2)-F{{l + x,)/2) dx^^dx^. (10) 

Let 

G(xi, x„) = F(x„) - Fix,) + F (x„/2) - F ((1 + xi)/2) . 

The integral in Equation (fTO|) is critical at (a;i,a;n) = (0, 1), since G(0, 1) = 1 and for (a;i,x„) e (0, 1)^ the 
integral converges to as n — > oo. So we make the change of variables zi = xi and z„ = 1 — x„, then 
G(xi,a;„) becomes 

G(zi, 2„) = i^(l - z„) - F(zi) + F ((1 - z„)/2) - F ((1 + ^i)/2) , 
and Equation (flO)) becomes 

Pn{F)K [ [ ?i(?i-l)/(zi)/(l-z„)[G(zi,z„)]""' dz„dzi. (11) 
Jo Jo 

The new integral is critical at (2i,z„) = (0,0). Note that — g^^ dz '= " ^ *-* ^"-"^ ^^-'^ i'tS > 1. Let := 
^^^gfe^l(o+,o+) =/«(0+)+2-(^+i)/« and;9, ""g:i""^ l(o^o.) = /(^■)(l-)+2-0-+i) /O") (i"). 

Then by the hypothesis of the theorem, we have = and (i^) = for alH = 0, 1, . . . , (fc — 1); and 

/3j = and ^ = for all j = 0, 1, ...,(£ — 1). So the Taylor series expansions of f{zi) around 

zi = 0^ up to order k and /(I — z„) around z„ = 0^ up to order £, and G{zi,Zn) around (0^,0^) up to 
order (fc + 1) and {£ + 1) in zi, z„, respectively, so that (zi, z„) £ (0, e)^, are as follows. 



/(^i) = ^/'''(o+) -'1 + o (.r^) ; /(I - .„) = 4 + o (4+^) ; 



ctfe fe+1 . ( — 1) 01 e+1 
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Then substituting these expansions in Equation (|lip . we obtain 



PniF) 



n{n — 1) 



JO 



kl 



/W(0+)zi^ + O(zf+^) 



4 + 0(4+^ 



«fc fc+i _ (-1)^/3^ e+i 
(fc + 1)! 1 (£ + 1)! 



71-2 



dzndzi. 



Now we let zi ~ wn i/C^'+i). z„ = f n 1/(^+1), and ly = min(fc,£) to obtain 



Pn{F) 



1 / Q!fc 



n (n—l) 



n (n — 1) 



(^+1)! 
(-1)^ 



-1 n-2 

/(fc)(0+)/W(l-)y;fc/ + O („-(2fc+3)/(fe+i)^ + o (n-(2«+3)/(^+i) 



O ('„-2(fe+2)(£+2)/((fe+l)(.+l))\ ^ _ 1 r^fe (1)1^ ^i+ll Q / _(.+2)/(.+ l)\ 

V y nl{k + iy. (£+1)! J V J 

letting n ^ oo, 

CO poo ( -. \f 

exp 



dvdw. 



/o Jo 



(fc + 1)! (£ + 1)! 

/(fc)(0+)/W(l-)(-l)^(fc + l)!(£ + l)! _ /(^)(0+)/(^)(l-) 



dvdw 



k\e\ {-lYik + !){£ + 












/W(0+)/W(l-) 




/('^)(0+) + 2 


-(fe+i) j-(fe) 1^1+^ 




/w(i-)+2-(^+i)/w (i y 



(12) 



as n ^ oo at rate 0{c{f) ■ n where c(/) is a constant depending on /. 



For the general case of 3^ = {vii y2}, the transformation (f){x) 



y2-yi 



maps (yi, y2) to (0, 1) and the trans- 



formed random variables <j>{Xi) are distributed with density g{x) — (y2 — yi) / (Oj !)■ Substituting 

f{x) by 5(2;) in Equation (|12|) . the desired result follows. ■ 

Note that 

• if min(/('=)(y+),/W(y-)) = and min (z^'^') , /^^^ then ^ as 

n 00, at rate 0{c{f) ■ n^™) and 



if min(/(^-)(y+), f^'\y^)) ^ and Z^'^) (^^^^) = /^'^ (^^^ ) = then p„(F) 
at rate 0(c(/) • n"™). 



1 as n ^ 00, 



For example, with F — W(yi,y2), in Theorem 15.11 we have k = £ = 0, f{yi) = /(yj ) = / 



(yi+y2)+^ 
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+y2)" 



= l/(y2 — yi)- Then liijin^oo Pn{F) — 4/9, which agrees with the result given in Equation 

Example 5.2. For F with density f{x) ^ {x + 1/2) l(0 < a: < l), wc have fc = £ = 0, /(0+) = 1/2, 
/(1-) = 3/2 and / = / (^i") = 1. Thus hm„^oo = 3/8 = 0.375. The numerically computed 

(by numerical integration) value of Pn{F) with n = 1000 is piooo(-P) ~ 0.3753. □ 

Remark 5.3. Let pp := Imin^ao Pn{F)- Then the finite sample mean and variance of 'jniF) are given by 
1 +pn{F) and Pn{F) (1 —pn{F)), respectively; and the asymptotic mean and variance of 'jniF) are given by 
1+Pf and pp (1 — pp), respectively. □ 

Remark 5.4. In Theorem 15. H we assume that and /W(-) are bounded on (yi,y2). Suppose either 

or f^^\-) or both are not bounded on (yi,y2) for fc,Z > 0, in particular at yi,(yi +y2)/2,y2, for 
example, lim^^y+ f'^^\x) = oo. Then we find p{F) as 



p{F) 



lim 

(5^0+ 







[/(fc)(yi+<5) + 2-(^+i) ((^^7^' +(5) 







□ 



Example 5.5. Consider the distribution with density function f{x) 



1(0 < a: < 1). Note that 



3^2 = {0, 1} and f{x) is unbounded at x G {0, 1}. See Figure [1] (left) for the plot of f{x). Instead of /(x), we 
consider g{x) = 2 arcsin'(l-2i5) -"-(^ < X < 1 — 5) with cdf G{x). For g{x), we have fc = £ = in Theorem 16.31 
and then lim„_>oo Pn (^") = limj^o+ linin^oo Pn(G) = 1 using Remark 15.41 The numerically computed value 
of piooo(i^) is Piooo(^^) ~ 1.000. □ 




Figure 1: Graph of the density in Examples 15.51 (left) and 15.71 (right). 
Remark 5.6. The rate of convergence in Theorem 15.11 depends on /. From the proof of Theorem 15.11 it 
follows that for sufficiently large n, 



Pn{F) 



/W(y+)/W(y-) 



c{f) 



f^^Hyt) + 2-(fc+i) /W {^-^^^^ .f'^^Ky2) + 2-(^+i) /W (^MiM j 



where 



c(/) 



„ c'^TT r i_ c ^TTi y ( k+2\ 



(fc + l)(£+l)s 
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with T{x) = cxp(-t)t(^-i) dt and 

- = ;;4r /""'(y^) /'"(y.-), ^ (/"'(y.-) + 2-.'«. /» (M") ) , 

provided the derivatives exist. □ 

Example 5.7. Consider the distribution with absolute sine density f{x) = tt/2 \ sin(2 7r2;)| 1(0 < x < 1). See 
Figured (right) for the plot of /(x). Then 3^2 = {0, 1} and since /(0+) = / (^i^) = and /(l") = / (^i") = 

and /'(O) = /' (i^) = TT^ and /'(l") = /' (^i^^ = we apply Theorem O with k ^ £ = 1. Then 

lim„^ooPn(^) = 16/25 = 0.64. The numerically computed value (by numerical integration) of piooo(^) is 
Piooo(J^) « 0.6400. □ 



The distribution of jn{F) depends on the distribution of r{Xi) = min{d(Xi,yi), c?(Xi,y2)). Based on 
this, we have the following symmetry result. 

Proposition 5.8. Let Fi and F2 he two distributions with support S{Fj) C (71,72) for j = 1,2 such that 
Fi{yi + x) = 1 — F2 (y2 — x) for all x E (0, y2 — yi) (hence /i(yi + x) ~ f2 (y2 — x) ). Also, let he a set of 
iid random variables from Fj for j = 1,2. Then the distributions of ^n{Fj) are identical for j = 1,2. 



Proof: By the change of variable X = (p{U) = y2 — yi — ?7 for t/ € (0,y2 — yi), we get F2 (yi + u) = 
1 — Fi(y2 — u). Furthermore, <p(u) transforms ri{X^, Ny) into Ti{X^,Ny) for A"^, so P(7„(Fj) = 2) are 
same for both .? = 1, 2. Hence the desired result follows. ■ 



Below are asymptotic distributions of 7n(-F) for various families of distributions. Recall that pp — 
liuin^oo Pn{F) = linin^oo -P(7ri(-F') = 2). The asymptotic distribution of jn{F) is 1 + B ernouUif j PF) . For 
the piecewise constant functions in Section 14.2.11 Theorem 15 . 1 1 applies . See Section 6.1 in CevhanI (|2004 ). 

Example 5.9. Consider the distribution F with density /(•) which is of the form 

f{x) = {ax + b) I(.T G (0, 1)) with |a| < 2, 6 = 1 - a/2. 

So fc = £ = and /(0+) = 6, /(l") = a + 6 and / = / (^i") = a/2 + 6. Then by TheoremO we have 

4 - a^ 

lim Pn{F) = 2 ~- PF{a)- 

Note that pf{ci) € [0,4/9] is continuous in a and decreases as \a\ increases. If a = 0, then F = U{0, 1), and 
PF{a = 0) = 4/9. Moreover, pria = ±2) = 0; that is, for a = ±2, the asymptotic distribution of jn{F) is 
degenerate. □ 

Example 5.10. Consider the normal distribution A/'(/i, cr^) restricted to the interval (0, 1) with /i e M and 
(T > 0. Then the corresponding density function is given by 

f{x,fi,a) = K (^^) cxp \ 1(0 < a; < 1), 
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where k = [$ {^-^) — $ ("7^)] ^ ^it^i '^'l') being the cdf of the standard normal distribution A/'(0, 1). Note 
that fc = £ = 0, then by Theorem [Q 



hm PniF) 



(2 + exp(^))(2 + exp(y#)) 



Observe that f) G [0, 4/9) is continuous in /x and a and increases as ct increases for fixed ^. Furthermore, 

for fixed /i, limo-^oo cr) = 4/9 and Umg-^o o") = 0. For fixed a > 0, lim^^^oo o") = 0, 

Pf{p-i c) decreases as \fi — 1/2| increases, and ppifJ-, f) is maximized at /i = 1/2. □ 

Example 5.11. Consider the distribution F with density /(•) which is of the form 

f{x) = 2%q + l) [a;-? 1(0 < X < 1/2) + (.T- l/2)« 1(1/2 < X < 1)] with g e [0,oo]. 

See Figure [2] (left) with q = 2. Since /(0+) = / = 0> "^c can apply Theorem 15. II with k = q and / = 0. 

Then /(«)(0+) = {q + 1)! 2?, /(l") = (<? + 1), / (5") = (<? + 1), and = (g + 1)! 2«. By Theorem 

29+2 



5.11 we have 



lim PniF) 



— =:pF{q)- 



3(1 + 29+1) 

Note that priq) G [4/9, 2/3] is a continuous increasing function of q. li q = 0, then F = U{Q, 1). □ 
Example 5.12. Consider the distribution F with density /(•) which is of the form 

fix) = ((5 + 12(1 - S)x^) 1(0 < a: < 1/2) + {6+12(1-6) {x-l/2)^) l(l/2 < x < l) with 6 e [0,1]. 

See Figure[2]with 6^0 (left) and 6 = 2/3 (right). Since /(0+) = 6, f (1") = {3-2 6), f (^^) = (3-2 6) 



and / 



1 + 



= 6, for 6 G (0, 1]} we have fc = £ = and so by Theorem 15.1 

lim PniF) = 4/9 for 6 G (0,1]. 



Note that if 5 = 1, then F = UiO, 1). For 6 = 0, we can apply Theorem 1 5 . 1 1 wit h k = 2 and I = 0. Hence we 
get ppiS = 0) = 16/27. Observe that in this example, jniF) has two distinct non-degenerate distributions 
at different values of (5. □ 




Figure 2: Left plot is for the density in Example 15.111 with q ~ 2 or for the density in Example 15.121 with 
6 = 0. Right plot is for the density in Example 15 . 1 21 with 6 = 2/3. 
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Remark 5.13. If, in ThcorcmO wc have f ''\0+) = (i"^) and f^'^Hl^) = f^'^^ (5 ), then 

lim Pn{F) : 



1 \ / 1 



1 + 2-^^+1) J \1 + 2-(^+i) 

In particular, if k = £ = 0, then hm„_>oo Pn (.F) = 4/9 (i.e., jniP) and 7„(Z^(0, 1)) have the same asymptotic 
distributions). □ 

Example 5.14. Beta(i/i, 1^2) with 1^1,1^2 > 1- The density function is 

/(z, ^1,^.2) - ''^'''{^"''f 1(0 < X < 1) where /3(^i,z.2) = ^P^^- 

Then Um„^ooPn(Beta(!/i, 1^2)) = at rate 0(n"('"i+'"2-2)^ ^ l^I- p^^^^^p^) denote the P[^n{F) = 2) for 
F = Beta(i'i, V2). The numericaUy computed values of Pn{vi, V2) for n = 1000 are piooo(4, 1) = Piooo(l) 4) « 
0.000005, piooo(4, 2) = piooo(2, 4) < 0.00001 and piooo(2, 2) w 0.000001. □ 

Here is an example with general support (yi,y2). 
Example 5.15. Consider the distribution F with density /(•) which is of the form f{x) = ax + b with b = 
^^(1 - a(yi - y?)/2) and \a\ < Using TheoremO we obtain pp = i^^lZll]lzt - « (Vi^Vs) - 

(0,1), then = l-a/2 and pi^(a) = In both cases, ppio.) is maximized for the uniform case; i.e.. when 

a = 0, then we have ppia = 0) = 4/9. Furthermore, '^n{F) is degenerate in the limit when a = , 
since Pn{F) as n cx) at rate 0(n~-'^). □ 



For more detail and examples, see Sections 6.4 and 7.1 in lCevhanI (j20041) 



6 The Distribution of the Domination Number of „^-digraphs 

In this section, wc attempt the more challenging case of m > 2. For c < d in R, define the family of 
distributions 

^(M) {FxY ■■ {X^.Y,) - Fxy with support S{Fxy) {c,df C M^^ X, - Fx and - Fy] . 

We provide the exact distribution of ^{Dn.m) for J^(M)-random digraphs in the following theorem. Let 
[m] := {0, 1, . . . ,m - 1} and Of := {(ui, ...Ub) : = a : Ui& S, Vi}. Let = {yL,>2, • ■ • ,i^m} 

whose order statistics are denoted as Yi^j^ for j = 1, 2, . . . , m. Note that the order statistics are distinct a.s. 
provided Y has a continuous distribution. Let j{D^) be the domination number of the digraph induced 
by and (see Section |4]). Given Y(j) = y^) for j = l,...,m, let be the (conditional) marginal 
distribution of X restricted to X, — (y(j-i) , y(j)) for j ~ 1, . . . , (m + 1), n be the vector of numbers of X 
points Uj falling into intervals Tj. Let fyif) be the joint distribution of the order statistics of Xn> i-e., 
/^(y) = rijli /(yi) I(c < yi < . . . < ym < d), and fj.kiVj^yk) be the joint distribution of l?,i, l?fcv Then 
we have the following theorem which is a generalization of the main result of IPriebe et al. I (I2OOII) . 



Theorem 6.1. Let D he an J^/f{W) -random &n.m-digraph. Then the probability mass function of the domi- 
nation number of D is given by 

P{l{Dn,rn) = k) = j ^ ^ P(iV = n) C(fci , JT-i ) ((fc^+i , ) J]^ //(fcj , ) (y) dyi . . . dy„ 
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where P(N = ft) is the joint probability of nj points falling into intervals Ij for j = 1,2, ... , (m + 1), 
fcj e {0,1,2}, and 

C{kj, Tij) = max(l(nj = kj = 0), > kj = 1)) for j = 1, (to + 1), and 
77(fc„n,) = max(l(n, = k, = 0),I(n, > k, > 1)) • (^^O^f^-^) _ ^^^^^^.^^i(/c,=i) 
for J = 2, . . . , TO, and t/ie region of integration is given by 
^ {(yi,y2, ■ • ■ ,ym) e (c, d)^ : c < yi < y2 < . . . < y,„ < d}. 



Proof: For j{Dn.m) ~ li^'') ~ k, we must have j{D^) = kj for j = 1, . . . , (to + 1) so that 

fcj = k and X^jS^ "-i = definition, ©I'^j^^+i) is the collection of such n and since kj G {0, 1, 2} 

for all j = 1, . . . , (to + 1), 0jj is the collection of such k. We treat the end intervals, Xi and Tm+i, 

separately. The indicator functions in the statement of the theorem guarantees that the pairs Uj , kj are 
compatible for j G {1, (m + 1)}; that is, incompatible pairs such as (uj ~ 0,kj > 0) are eliminated. The C 
terms equal unity if (uj, kj) are compatible. Therefore we have 

^ ?n+l 
P{l{Dn,nr) = k) = ^ ^ P(7V = Tl) \{ Tj{k, , n,) ff{y) dyi . . . dy,n 

P m 

= P{N = ri) Y[ vikj,nj)Y[T]{kj,nj)f^{y)dyi...dy,n 

n,(m + l) k,(m + l) 

» m 

= Y Y PiN = n)C{ki,ni)Cikrn+i,n„i+i)Y[vikj,nj)f^{y)dyi...dyrn 



i=2 



where we have used the conditional pairwisc independence of ^{D^). The r/ terms arc based on the compat- 
ibility of pairs {uj, kj) for j = 1, . . . , (to + 1) and (Fj) = P{j{D^) = 2). ■ 

For n,m < oo, the expected value of domination number is 

m n 

E[7(^„.™)] = P (X(i) < 1^(1)) + P (^(„) > V)) + E E ^) E[7(i?^)] (13) 

where 

P{Nj =k) = 

/ /j-ij (yo-i)'yo)) -^^ (VO)) - (yo-i)) l - (-f^^ (VO)) - (yo-i))) dyo)dyo_i) 

c "'¥(3-1) 

andE[7(I?J")] = l+pfc(F,). 

Corollary 6.2. For Fx, y G ^(M) twii/i support S (Fx) H S{Fy) of positive measure, lim„_,oo E[7(D„^„)] = 

00. 



Proof: Consider the intersection of the supports S{Fx) H S{Fy) that has positive (Lebesgue) measure. 
For S{Y) \ S{X); i.e., in the intervals Ij falhng outside the intersection S{Fx) n S{Fy), the domination 
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number of the component is j{D^) ~ w.p. 1 but inside the intersection, j{D^) > w.p. 1 for infinitely 
many j. That is, 

n n 

E[7(^n.n)] = P < y(i)) + P(X(„) >!-(„)) + E E ^(^J- = ^) (^.)] 

j=2 fc=l 

n n n n 

> E E ^(^^- = ^) E[7JV, (F,)] = E E ^(^^- = (1 + (Fj)) 

j=2 k=l j=2 k=l 

> EEm = fc)>Em>i) 

j=2 k=l j=2 

« n (for sufficiently large n) 

where E[7Ar (-F})] = (1 +pn- (Fj)) follows from the fact that (Fj) ~ l + Bernoulli(p7v {Fj)) from Theorem 
14.21 Furthermore, P{Nj > 1) ~ 1 for sufficiently large n. Then the desired result follows. ■ 

Theorem 6.3. Let Dn_m be an J^(M.) -random S^n m-digraph. Then (i) for fixed n < oo, limm^oo l{Dn,m) = 

n a.s. (a) for fixed m < oo, lim,i_>oo l{Dn,m) = m + 1 + X^Jli Pjj where Bj ~ Bernoulli (pi^^ ) where = 
stands for equality in distribution. 



Proof: Part (i) is trivial. As for part (ii), first note that Nj ^ oo as n — > oo for all j a.s., hence 
lim„^oo l{D^) = lim,i_»oo 'y{D™'^^) = 1 a.s. and lim„^oo l{D'') = 1 + Bernoulli(pi7v ) a.s. for j = 2, . . . , to 
where 



PF,^ / H* (yo-i),y(j)) fj-ij (yo-i)-y(j)) dyu)dy (j-i) 

with H* (y(j-i) 1 y(j")) = liiirnj^oo(Pnj (-Fj)) which is given in Theorem 16.11 for Fj with density fj whose 
support is (y(j-i) , y(j)) • Then the desired result follows. ■ 

So far, ym is assumed to be a random sample, so P{"i{Dn.m) = k) includes the integration with respect 
to /y(y) which can be lifted by conditioning. Conditional on y,n = {y(i), • ■ • ,y(m)}j by Theorem 16.11 we 
have 

Ph{F>7i,r,i) ^ k) = ^ ^ P{N = n)C{ki,ni)C{kjn+i, n,n+i)Y[v{kj,n.j), 

where C,{kj,nj) and ri{kj,nj) are as in Theorem 16.11 and the expected domination number E[7(Z)„^m)] is 
as in Equation ^ with P{Nj ^ k) ^ [Fx (yo)) - Fx (yo-i))]'' [l - {Fx (vo)) - Fx (y(j-i)))] ""^ and 
hm„^oo 7(-D„,m) = m + 1 + X^Jli Pj^ where Bj ^ BcrnouUi(pF3 ) with pp^ hm„^.^oo Pn^ {Fj)- 

Let Fx be a distribution with support S{Fx) C (0,1) and density fx{x). Conditional on 3^,„ = 

{y(i) , • ■ • , y(m) } , let Fj be the distribution with density fj (.t) = ^ ^ ( y^~-yj'',J for j = 2, . . . , m, 

and S{Fj [x)) is non-empty for G {1, (to+ 1)}. By this construction, the independence of the distribution of 

jniFj) from 2j obtains; i.e., 7„(i^j) Jn{Fx) for all j G {1, . . . , (to + 1)}. Now consider the family J^(R) 
defined as 

JiiiiR) ~ {^^x.y : {X,,Y,) ^ ^^^.y, '^U{c,d) for (c, d) C M, and X,;|3^^ I'' F,}. 
Clearly J^(M) C 
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Corollary 6.4. Suppose Fx.y G ^(M). Then 



where C,{kj,nj) and rj{kj,nj) are as in Theorem \6.1 



Note that if in addition, Pf,{X e Ij) = G X,) for all j, then P{N ^ n) = ("+") \ since each 

^ ^ ®rr(m+i) occurs with probability ("^J") ^- Moreover, F =U{c, d) is a special case of Corollary 16.41 For 
n, m < (30, we have the explicit form of Pn^ (Fj) for Fj with piecewise constant density fj. 

Here are some examples which are generalized from piecewise-constant densities so that now the distribu- 
tion of j{D^) is independent from the support (y(j_i),y(j)). Hence Corollary 16.41 applies to these examples. 

Example 6.5. Let Uj :~ i'^'-'-'->^^'-^>) g^j^j y^^.^ _ y(j-i)- 

• If /(•) is of the form 

then Pn{F) is as in Equation ([7]). 

• If /(•) is of the form 

^^^^ = (1 -^(^ ^ (yO-i)'"j - '^^j) U K + (5wj,y(j))) with^e [0,1/3], 

then Pn{F) is as in Equation ([8]). 

• If /(•) is of the form 

/(a,) = ii±^ I (a; e (y(^._,),u^.)) + ii_^l(a;e h,y(,))), 

then Pn{F) is as in Equation 

• If /(•) is of the form 

f{x) ^ h{x)l{x £ (y(j_i),ij)) +/2(a;)I(a;e [t^, w^)) + /g (a:) I (a; S K-,yo))) 

where t, = ^^iUli^, = ^Jl,2±n^^ /^(^) ^ M , (^.) ^ and /a (x) = (i±^, then p„(F) 

is as in Example 14.91 □ 

Theorem 6.6. Let D he an jy/^{M.) -random !3n.m-digraph with the additional assumption that Pp.{X G 
Xj ) = Pu {X e Xj ) for all j . Then 

„r . „ M 2 n n! m (m — 1) (n + m — i — 1)\ , , , 

n-\-m n + ml! ^ — ^ n — z ! 

where p,{F) ^ P{j{D^^2) = 2). 
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Proof: Similar to the Proof of Theorem 4 in iPriebe et ahl (j200ll ) 



Furthermore, from CoroUarv 16.21 we have E[7(Z?„.„)] oo as n — s- oo. 

Theorem 6.7. Let Dn^m be an Ji^{M.)-random &n,m-digraph. Then (i) for fixed n < oo, hm,„_>oo7(^ri 
n U.S. (ii) for fixed m < oo, hm„_ 
PF = hm„_»oo P{l{D„,2) = 2). 



• l{Dn,m) = m + 1 + B , where B ~ Binomial{m — 1,pf) where 



Proof: Similar to the Proof of Theorem 5 in iPriebe et al.l (j200l[ ) 



Remark 6.8. Extension to Multi-dimensional Case: The existence of ordering of points in R is crucial 
in our calculations. The order statistics of partition the support (c, d) into disjoint intervals a.s. which 
can also be viewed as the Delaunay tessellation of R based on 3^,„. This nice structure in M avails a minimum 
dominating set and hence the domination number, both in polynomial time. Furthermore, the Fi-region 
is readily available by the order statistics of Xn', also the components of the digraph restricted to intervals 



for X, y in distinct intervals, 
does not have a nice ordering 



Ij (see Section HJ are not connected to each other, since Ny{x) n Ny{y) ~ 
The straightforward extension to multiple dimensions (i.e., M'' with d > 1 
structure; and 3^,„ does not readily partition the support, but we can use the Delaunay tessellation based 
on y,n- Furthermore, in multiple dimensions finding a minimum dominating set is an NP-hard problem; 
and Fi-regions are not readily available (in fact for nj > 3, complexity of finding the Fi-rcgions is an open 
problem). In addition, in multiple dimensions the components of the digraph restricted to Delaunay cells 
are not necessarily disconnected from each other, since Ny{x) n Ny{y) ^ might hold for x, y in distinct 
Delaunay c ells. These have mot i vated us to generalize the proximity map Ny in order to avoid the difficulties 
above. See lCevhan and Priebd (j2003l . 120051 ). where two new families of proximity maps are introduced, and 
the generalization of CCCD are called proximity catch digraphs. The distribution of the domination number 
of these proximity maps is still a topic of ongoing research. □ 



7 Discussion 



This article generalizes the main result of IPriebe et al.l (j200lh in several directions. iPriebe et all (|200ll ) 



provided the exact (finite sample) distribution of the class cover catch digraphs (CCCDs) based on 
and y^a both of which were sets of iid random variables from a uniform distribution on (c, d) C R with 
— cx) < c < d < oo and the proximity map i^y{x) := B{x,r{x)) where r{x) := miny^y^ d{x,y). First, 
given 3^2 = {yi,y2} C R, we lift the uniformity assumption of Xn by assuming it to be from a non-uniform 
distribution F with support S{F) C (yi,y2). The exact distribution of the domination number of the 
associated CCCD, 7„(F), is calculated for F that has piecewise constant density / on (yi,y2). For more 
general F, the exact distribution is not analytically available in simple closed form, so we compute it by 
numerical integration. However, the asymptotic distribution of 7n(-F) is tractable, which is the one of the 
main results of this article. Unfortunately, the distribution of jn{F) depends on 3^2, hence the distribution of 
the domination number of a CCCD, j{Dnjn)j for <Y„ and ym with m > 2, for general F includes integration 
with respect to order statistics of X„. We provide the conditions that make j(Dn,m) independent of y,n- 
As another generalization direction, we also devise proximity maps depending on F that will yield the 



distrib ution identical to that of 7Ti(i^(yi, y2))- Our set-up is more general than the one given in IPriebe et al 



(j200ll ). The definition of the proximity map is generalized to any probability space and is only assumed to 



have a regional relationship to determine the inclusion of a point in the proximity region. 

The exact (finite sample) distribution of jniF) characterizes up to a special type of symmetry (see 
Proposition 15. 8p . Furthermore, this article will form the foundation of the generalizations and calculations 
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for uniform and non-imiform cases in multiple dimensions. As in ICevhan and Priebd I 2005f) . we can use 
the domination number in testing one-dimensional spatial point patterns and our results will help make the 
power comparisons possible for large families of distributions. 
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